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ABSTRACT 

The incidence of strong Mgll systems in gamma-ray burst (GRB) spectra is a few times higher than 
in quasar (QSO) spectra. We investigate several possible explanations for this effect, including: dust 
obscuration bias, clustering of the absorbers, different beam sizes of the sources, multiband magnifica- 
tion bias of GRBs, association of the absorbers with the GRB event or the circumburst environment. 
We find that: i) the incidence rate of Mgll systems in QSO spectra could be underestimated by a 
factor 1.3-2 due to dust obscuration; ii) the equivalent- width distribution of the Mgll absorbers along 
GRBs is consistent with that observed along QSOs thus suggesting that the absorbers are more ex- 
tended than the beam sizes of the sources; Hi) on average, GRB afterglows showing more than one 
Mgll system are a factor of 1.7 brighter than the others, suggesting a lensing origin of the observed 
discrepancy; iv) gravitational lensing (in different forms, from galaxy lensing to microlensing) can bias 
high the counts of Mgll systems along GRBs if the luminosity functions of the prompt gamma-ray 
emission and of the optical afterglows have a mean faint-end slope approaching -5/3 - -2; v) some of 
the absorbers can be associated with the circumburst environment or produced by supernova rem- 
nants unrelated to the GRB event itself but lying in the same star-forming region. With the possible 
exception of magnification bias, it is unlikely that one of these effects on its own can fully account 
for the observed counts. However, the combined action of some of them can substantially reduce the 
statistical significance of the discrepancy. 

Subject headings: quasars: absorption lines - gamma-rays: bursts 
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1. INTRODUCTION 

Magnesium is an a-process element produced by red- 
giant stars and dispersed in the interstellar medium by 
supernova explosions and stellar winds. In the red- 
shift interval 0.3< z < 2.2, the Mgll doublet (2796, 2804 
A) produces absorption lines in the optical spectrum 
of background sources. The strong absorption systems 
(with equivalent width W > 0.3 A) are thus believed 
to be good tracers of metal-enriched gas associated with 
galaxies (e.g. Bergeron & Boisse 1991; Steidel, Dickinson 
& Persson 1994). Typically these systems have multiple 
velocity components (Churchill and Vogt 2001). 

Recent high-resolution imaging studies of quasar fields 
provided evidence that strong Mgll absorption is pro- 
duced in patchy gaseous envelopes, up to impact param- 
eters of 80 hr 1 kpc, surrounding galaxies of different mor- 
phological types (Churchill, Kacprzak & Steidel 2005 and 
references therein) . The completion of large and homo- 
geneous quasar (QSO) samples as the Sloan Digital Sky 
Survey (SDSS; York et al. 2000) and the Two-degree 
field quasar survey (2Qz; Boyle et al. 2000; Croom et al. 
2004) allowed accurate statistical studies. The amplitude 
of the cross-correlation function between Mgll systems 
and luminous red galaxies suggests that absorbers with 
W > l A are hosted within dark-matter halos with 
characteristic masses of 10 11_12 M Q and could be as- 
sociated to galactic superwinds (Bouche et al. 2006). 
Using the SDSS Early Data Release, Nestor, Turnshek 
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& Rao (2005) identified over 1,300 Mgll doublets with 
W > 0.3 Aand measured their equivalent width distribu- 
tion over the redshift range 0.366 < z < 2.269. Similarly, 
Prochter, Prochaska & Buries (2006a) found nearly 7,000 
Mgll systems with W > 1 A in the spectra of 50,000 
QSOs from the SDSS DR4. This corresponds to a red- 
shift path density dN/dz ~ 0.24 at z = 1. These results 
have been extended to lower redshift (and weaker col- 
umn densities) by Nestor, Turnshek & Rao (2006). This 
study suggests that the gas clouds associated with weak 
Mgll systems are a physically distinct population from 
those producing strong Mgll absorbers (see also Nestor 
et al. 2005). 

Gamma-ray bursts (GRBs) with bright optical after- 
glows can also be used as background sources. Somewhat 
surprisingly, Prochter et al. (2006b, hereafter P06) iden- 
tified 14 Mgll systems with W > 1 A along 14 GRB lines 
of sight (for a total redshift path of 15.5 at a mean red- 
shift z = 1.1). This corresponds to dN/dz = 0.90^;^ 
(symmetrical 99% Poisson 4 confidence interval), an in- 
cidence rate that is a few times larger than that in- 
ferred from the QSOs of the SDSS data set. A simi- 
lar (but much less statistically significant) discrepancy, 
dN/dz = 0.62^49 (with an additional - 30% uncer- 
tainty due to the redshift path), has been reported by 
Stocke & Rector (1997) in the optical spectra of 21 radio- 
selected BL Lacertae objects. However, recent redshift 
determinations (see e.g. Sbarufatti, Treves & Falomo 
2005) showed that 1 out of the 5 strong systems found 
by Stocke & Rector (1997) is associated with the host 
galaxy. This reduces the cosmological incidence rate to 

4 Clustering of the absorbers produces super-Poisson fluctua- 
tions, this effect will be discussed in i(2] 
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dN/dz = 0.49±5:2S- 

In this work, we critically review some possible ex- 
planations of the discrepancy between the counts of 
Mgll systems along GRBs and QSOs. The paper is or- 
ganized as follows. In fj2[ we estimate the importance of 
dust-obscuration bias along lines of sight towards QSOs. 
The effect of different beam sizes between GRBs and 
QSOs and the solution proposed by Frank et al. (2006, 
where QSO beams are assumed to be larger than GRB 
ones by a factor of 2) are discussed in <JJJ The roles of 
gravitational lcnsing and magnification bias of GRB af- 
terglows are presented in $4] Finally, in we discuss the 
possibility that some of the Mgll systems along GRBs 
are produced in the circumburst environment. Our re- 
sults are summarized in £|6] 

2. DUST OBSCURATION BIAS 

The presence of dusty absorbers along the line of sight 
could obscure the optical light from background QSOs 
and produce a selection bias in magnitude-limited sam- 
ples (Ostriker & Heisler 1984; Heisler & Ostriker 1988; 
Fall & Pei 1993). If the obscuration bias is important, ra- 
dio selected QSOs (which are unaffected by dust) should 
present a larger number of absorbers on average. Re- 
cent studies did not find strong evidence for dust red- 
dening and extinction (Ellison et al. 2001; Akerman et 
al. 2005; Jorgenson et al. 2006). However, radio samples 
are very small and cannot lead to definitive conclusions. 
Even though the number counts of absorbers in radio 
and optically selected QSOs are in agreement, la sta- 
tistical uncertainties are consistent with 60% of damped 
Ly-a systems being missed in optical magnitude-limited 
surveys (Ellison et al. 2004). 

Could dust obscuration bias of the QSO samples ex- 
plain the discrepancy with the incidence of Mgll systems 
in GRB spectra? A few dusty systems with color excess 
E(B — V) ~ 0.1 have indeed been found (Junkkarinen et 
al. 2004, Ellison et al. 2006). On the other hand, statisti- 
cal studies give contrasting evidence for dust obscuration. 
Characteristic values of E(B — V) = 0.06 — 0.1 have been 
reported for Znll and Call systems with large W (Vladilo 
k Peroux 2005; Wild et al. 2006). However, Murphy 
& Liske (2004) found no evidence for dust-reddening of 
QSOs by foreground damped Lyman-a systems. Simi- 
larly, York et al. (2006) could measure some appreciable 
reddening only in QSOs with very strong Mgll absorbers 
(W > 1.53 A). Moreover, radio and X-ray selected SDSS 
quasars do not appear to be more reddened than opti- 
cally selected ones. In both cases the typical color excess 
is E(B-V) ~ 0.01. 

All these results are puzzling and contradictory. In or- 
der to estimate the importance of dust obscuration in 
the SDSS sample of Mgll absorbers we have developed 
the following simplified model. We consider a popula- 
tion of QSOs at z = 2.3 with luminosities distributed 
according to the luminosity function determined by the 
SDSS (Richards et al. 2006). These QSOs will be ob- 
scured by foreground galaxies. To model the galaxy dis- 
tribution, we use 24 mock light cones extracted from the 
largest N-body simulation of the concordance cosmol- 
ogy performed so far, the Millennium run (Springel et 
al. 2005), in which galaxies are associated with dark- 
matter halos via semi-analytic modelling (Kitzbichler et 
al. 2006). The past light cone of an observer (i.e. an 
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Fig. 1. — The underlying number density of strong Mgll systems 
and the fraction of dusty absorbers that match the observed abun- 
dance of Mgll systems (solid) and the observed fraction of ab- 
sorbers along very reddened lines of sight (dashed) in the SDSS. 
The dotted lines mark the 0,1,2 and 3<r Poisson fluctuations for the 
observed density along GRBs. 



event with coordinates x, y, z, to) consists of all the paths 
of light that reach the observer at to- Basically, it con- 
tains all the galaxies that the observer could detect with 
an ideal experiment since they are connected by null 
geodesies with the observer himself. We shoot 10 3 ran- 
dom lines of sight for each light cone and count how 
many galaxies we find within a given impact parameter 
b = SOh- 1 kpc (Churchill et al. 2005). We then as- 
sociate a Mgll absorber to an intervening galaxy proba- 
bilistically with a covering factor / = 0.5 (Churchill et al. 
2005). We only consider galaxies with (unextincted) ab- 
solute magnitude Mb < M t h r and choose M t h r to match 
a given dN/dz. As a first case, we also assume that 
Mgll systems contain some dust with a bimodal distri- 
bution: a fraction /<j of them has E(B — V) = 0.1 (with 
a Small Magellanic Cloud (SMC) extinction curve) while 
all the rest has E{B — V) = 0.01. Dust obscuration can- 
not be discussed separately from gravitational lensing 
that boosts the luminosity of background sources. For 
this reason, we also compute the magnification due to 
all intervening mass concentrations along a given line of 
sight assuming that they are singular isothermal spheres 
(the resulting magnifications have a mean of 1.03 and an 
r.m.s. value of 0.04 but the distribution is very positively 
skewed, nearly 0.8% of the sources are amplified by more 
than 20%). 

Our results are summarized in Figure [T] The ab- 
scissa shows the true number density of the underlying 
Mgll systems (those revealed along GRB lines of sight), 
while the ordinate gives fd- the fraction of dusty ab- 
sorbers. The solid line indicates the parameter pairs for 
which SDSS would observe dN/dz — 0.24. Obscuration 
bias can decrease the observed density of absorbers by a 
maximum factor of ~ 3. In order to match the observed 
density of strong Mgll systems in SDSS, the number 
density of underlying absorbers cannot be higher than 
dN/dz ~ 0.7 but, in this case, all Mgll system should 
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be very dusty. This would be in contrast with obser- 
vations. York et al. (2006) found that only 14% of 
the Mgll systems with W > 0.3 A in the SDSS cata- 
log (and 36% of those with W > 2.4 A) lie along very 
reddened lines of sight (A(g — i) > 0.2). If this fraction 
amounts to ~ 20% for W > 1 A, the only way to rec- 
oncile our results with the observational data is the case 
where dN/dz ~ 0.33 (corresponding to Mthr — —21.65) 
5 and /d — 0.33. In this case, SDSS would miss 16% of 
the QSOs which are intrinsically brighter than its magni- 
tude limit (the effect of color selection is less important) 
and include 1% of the total number due to magnification 
bias (plus dust obscuration). Note that the median and 
mean halo masses of the Mgll absorbers are 3.0 x 10 11 
Mq and 7.5 x 10 11 M Q , in good agreement with Bouche 
et al. (2006). We explicitly checked that our results do 
not change substantially by using a Milky Way (MW) 
extinction curve for the dusty absorbers. 

In order to test if our results depend on the assumed 
distribution of reddening, we repeated our Monte Carlo 
simulations assuming a two-parameter Weibull distribu- 
tion for E(B — V) (for an SMC extinction curve). The 
cumulative Weibull distribution for a variable x > 
is C(x) = 1 — exp[— (ir/A) 7 ] and it is fully determined 
by the shape parameter 7 > and the scale parame- 
ter A > 0. This distribution is interesting because it 
can attain many different shapes based on the value of 
7. In particular, for 7 < 1 the corresponding proba- 
bility density function decreases monotonically and is 
convex; for 7 = 1 it becomes the exponential distribu- 
tion; for 7 > 1 it vanishes at x — and admits a mode 
at x = A(l— I/7) 1 / 7 ; for 7 = 2 it gives the Rayleigh 
distribution; for 7 < 2.6 it is positively skewed; for 
2.6 < 7 < 3.7 it closely approximates a normal distri- 
bution; for 7 > 3.7 it is negatively skewed. We find that 
only monotonically decreasing distributions for E(B — V) 
with 7 ~ 0.3 — 0.5 are compatible with the data. 

In this case, there are many solutions that match both 
the density of Mgll absorbers and the fraction of ab- 
sorbers along very reddened QSOs observed by SDSS (see 
Figure [2]). However, if we also require that the median 
observed color excess is of the order of 0.01 (as found in 
the SDSS by York et al. 2006) we have to assume that the 
underlying incidence rate is dN/dz ~ 0.30 (correspond- 
ing to M t hr — —21.7). This increases up to dN/dz ~ 0.45 
(Mthr ^ -21.4) if the 2Qz E(B-V) value (0.04, Outram 
et al. 2001) is adopted. 

In summary, we found that estimates of dN/dz based 
on QSO spectra are likely to be underestimated by a 
factor ~ 1.3—2 if dust-obscuration bias is important. 
If dN/dz ~ 0.35 (0.45), we find that there is a 0.15% 
(1.3%) chance to find 14 absorbers or more in a red- 
shift path of 15.5 because of random fluctuations. This 
shows that it is unlikely that dust obscuration bias fully 
explains the difference in the number of Mgll absorbers 
along GRBs and QSOs. Note that in our simulations, 
the scatter of the number of absorbers per line of sight is 
generally larger than expected for a Poisson distribution 
because the host-galaxies of the absorbers are clustered. 
If, on average, there are N absorbers per line of sight, 

5 Note that reducing either the maximum impact parameter b 
or the covering factor / would make the absorbers' host galaxies 
fainter for a given dN/dz. 
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Fig. 2. — Parameters of the reddening distribution of strong 
Mgll systems. Thick and thin curves respectively mark the set 
of parameters for which the observed abundance of absorbers and 
the fraction of very reddened lines of sight match the SDSS values. 
Solid, long-dashed, short-dashed and dotted lines refer to dN/dz = 
0.90,0.66,0.49,0.35 (i.e. to 0,1,2 and 3<r Poissonian fluctuations 
in the counts of GRB absorbers) respectively. The shaded area 
is bounded by the loci where the median observed color excess 
A[ln(2)] 1 / 7 is 0.01 (lower boundary) and 0.04 (upper boundary). 
Models for which thin and thick curves of the same type cross 
within the shaded area are consistent with all the observational 
constraints. 



the variance of the counts is a 2 = N (1 + N where £ 
is the mean two-point correlation function computed by 
averaging over pairs of points both lying within a narrow 
cylinder with radius b parallel to the line of sight. For 
b = 80 h" 1 kpc and a redshift interval 0.3 < z < 2.2, 
the galaxies hosting the absorbers have typically mean 
correlations of some percent (slightly depending on the 
mean galaxy luminosity). Whenever N^t 1, supcrPois- 
son fluctuations are then appreciable. A number, X, of 
independent lines of sight have to be combined together 
to obtain a redshift path of 15.5. This increases both the 
mean counts and the scatter by a factor of X so that the 
total excess of the scatter with respect to Poisson remains 

1 + N {; as for single lines of sight. In consequence, the 
scatter in the number of Mgll systems is nearly Poisso- 
nian along QSOs and superPoissonian along GRBs. Fu- 
ture data with increased statistics, might then be used 
to estimate the clustering amplitude of GRB absorbers. 

3. STATISTICS OF ABSORPTION LINES AND BEAM SIZE 

Frank et al. (2006, herafter F06) proposed a geometri- 
cal explanation for the different incidence rate of strong 
Mgll systems in QSO and GRB lines of sight. They ar- 
gued that the difference in the statistics can be readily 
explained if the size of the QSO beam is approximately 

2 times larger than the GRB beam, provided the size 
of the Mgll absorbers is comparable to the GRB beam 
and of the order of ^ 10 16 cm 2 . This is in contrast with 
a number of theoretical and observational estimates sug- 
gesting that the typical QSO beam is a few times smaller 
than the GRB beam (see Section 4 in F06 and references 
therein). We show here that the solution worked out by 
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F06 did not include an additional effect which changes 
the outcome of their model. 

In what follows r is the comoving radial distance and 
D a is the angular diameter distance. For simplicity, let 
us consider a population of gas clouds with comoving 
number density n{z), proper cross-section a and equiv- 
alent width distribution /(Wo) (normalized to unity). 
The clouds can be detected as absorption lines in the 
spectrum of background continuum sources with an- 
gular beamsize fib- If the beamsize is much smaller 
than the solid angle subtended by the gas clouds, i.e. 
fib *C cr/D 2 (z), the number density of absorption lines 
per unit redshift and equivalent width is 



-j^—j^ = an(l + z) 2 — /(Wq) 



(1) 



On the other hand, if the clouds are much smaller than 
the beam of the background source, i.e. fib ^> cr/D^z), 
two new effects need to be accounted for. First, the ob- 
served equivalent width is proportional to the covering 
factor of the cloud (i.e. to the fraction of the beam area 
covered by a single absorber) W bs = 17 /[^hD 2 {z)\ Wo = 
[er/<7b(z)] Wo- Second, the mean number of absorbers 
scales proportionally to the area of the beam (note that 
this effect was apparently neglected by F06) so that: 



d 2 N 



dW ( 



— =a b n(l + z) - 

obs dz dz a 



f(^W ohs 
\ a 



(2) 



The last possibility is obtained when Ob — cr. In this 
case, the number of absorption lines is 

^c(^+V^)Ml + z) 2 ^, (3) 
dz dz 

and Monte Carlo simulations must be used to derive the 
corresponding distribution of W bs that depends on the 
relative positions and shape of beam and absorbers. 

Let us now consider two different classes of sources 
(for instance GRBs and QSOs), and denote their rela- 
tive beamsize by x — f2c/^Q- At every redshift, from 
equation ([2]) we derive 

d2NG : {w G ,z) = x 2 4^rAw Q = xW G ,z) . (4) 



dW G dz 



dWqdz 



Therefore, to find four times more absorbers with 
W > 1 A along GRBs than along QSOs, the relative 

beam size must satisfy the equation 4 ^WQdz ^Wq = 



f?L ffik d W G = -12 dW Q . Adopting the fit 

for the equivalent width distribution derived either by 
Nestor et al. (2005) or by Prochter et al. (2006a) it is 
easy to show that this equation does not admit solutions. 

This applies to small absorbers and only approximates 
the case in which absorbers and beams have similar 
sizes. We performed a number of Monte Carlo simu- 
lations varying the beam and absorber sizes. For GRB 
afterglows we both considered disk and ring geometries 
(with varying thickness) while we only considered disk- 
like QSO beams. In no case could we reproduce the ob- 
servational results. If the F06 solution holds, the distri- 
bution of Mgll equivalent widths along GRBs should be 
flatter than in QSOs (this is also evident in our Monte 
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Fig. 3. — The cumulative probability distribution of the 
Mgll equivalent widths along GRBs listed in P06 (solid histogram) 
is compared with the fit for QSO absorbers given by Prochter et 
al. (2006a). The two distributions are compatible at the 95.7 % 
confidence level in the Kolmogorov-Smirnov sense. 



Carlo simulations). However, we found that the equiva- 
lent widths listed in P06 are perfectly consistent with 
the distribution derived from SDSS QSOs (Figure [3]). 
The Kolmogorov-Smirnov test suggests that the observed 
equivalent widths for GRB absorbers are compatible (at 
the 95.7 % confidence level) with being a random sam- 
pling of the probability distribution derived from SDSS 
QSOs. 

There are some other observational results that the so- 
lution proposed by F06 cannot explain: i) along QSOs 
there are no apparently unsaturated Mgll absorption 
lines with a doublet ratio (1:1), as would be ex- 
pected if small saturated systems were being diluted 
by a larger QSO beam; ii) an apparently unsaturated 
Mgll absorption system with a doublet ratio (1:1) has 
been detected along GRB030226 (Shin et al. 2006) in 
the spectrum of GRB030226 thus suggesting partial cov- 
ering of the GRB beam; in) if the strong Mgll systems 
resemble intercloud medium of the MW their sizes are 
likely to be in the range 1-1000 pc (e.g. Kobayashi et al. 
2002; Rauch et al. 2002; Churchill et al. 2003; Ding et 
al. 2003) which are difficult to reconcile with the sizes 
proposed by F06. 

4. STATISTICS OF ABSORPTION LINES AND LENSING 

Only a few spectra of optical GRB afterglows have 
been taken. Is it possible that the sample is heavily bi- 
ased towards lines of sight intersecting a large number of 
absorbers? Gravitational lensing due to intervening ma- 
terial could in principle boost the afterglow luminosity 
by amplifying its beam size. This would make optical 
spectra easier to take. 

Combining the data by P06 with those by Nardini et 
al. (2006), we find that afterglows with more than one 
Mgll absorbers are, on average, a factor of 1.7 brighter 6 

6 In terms of their optical luminosity measured 12 hours (in the 
GRB rest frame) after the trigger. 
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than the others. According to a t-Student test, there 
is a 10% chance that this difference is due to random 
fluctuations. If this result is strengthened by increased 
statistics (currently only four afterglows showing more 
than one strong Mgll absorber are known), it is then 
likely that some form of lensing caused by mass concen- 
trations associated with the absorbers themselves is the 
cause of their increased detection probability. 

The galaxies hosting the absorbers represent obvious 
lens candidates. However, since Mgll absorbers are seen 
at large distances from galaxy cores, it is unlikely that 
they can produce large magnifications. Our Monte Carlo 
simulations presented in S}2] show that lines of sight with 
more than 2 Mgll absorbers have a mean magnification 
of 1.04 and only 2% of them are magnified by more than 
a factor of 1.2 (see also Menard 2005). Strong lensing is a 
very rare phenomenon. Only a few GRBs in every thou- 
sand detections are expected to be strongly lensed by 
galaxy-sized halos (Porciani & Madau 2001). Millilens- 
ing by 10 7 M Q haloes has an optical depth which is about 
1000 times smaller than the strong lensing one (Porciani 
& Madau 2000). 

Afterglow spectra are generally taken within the first 
few hours after the gamma-ray triggering event when 
most of the optical emission is expected to come from 
a narrow ring of radius r s ~ 4 x 10 15 (t/1 hr) 5 / 8 (l + 
zgrb) -5 / 8 cm which is expanding at super luminal speed 
on the sky (Waxman 1997). The beam size is thus 
comparable with the Einstein radius of compact so- 
lar mass objects at cosmological distances r E ~ 5 x 
1O 16 (M/M ) 1 / 2 and GRB afterglows can be efficiently 
microlensed by intervening stars and massive compact 
objects (MACHOs). When the source can be regarded 
as pointlike with respect to the lens, microlensing pro- 
duces a constant amplification depending on the im- 
pact parameter of the lens b (Loeb & Perna 1998). On 
timescales of days (observer frame), the magnification 
then increases and reaches a maximum value when the 
source size crosses the lens (r s = b). This rather sharp 
brightening is a characteristic signature of microlensing 
events. For large impact parameters, when the mean 
magnification is < 2, and for broader source rings this 
feature is less pronounced (see Figure 1 in Loeb & Perna 
1998) and could then be not easily detected observation- 
ally. For sources at z ~ 2, the microlensing optical depth 
is r ~ 0.65 fijjj where fi^i denotes the present-day density 
of microlenses in units of the critical density of the Uni- 
verse (Baltz & Hui 2005). Even in the most optimistic 
assumption that 20% of the dark matter is in MACHOs, 
only a few percent of the GRB afterglows should be af- 
fected by microlensing. In summary, unless the intrinsic 
luminosity function of GRBs (and of their afterglows) 
is extremely steep and magnification bias plays an im- 
portant role (see below), a lensing explanation of the in- 
creased Mgll absorption along GRB lines of sight should 
be regarded as very unlikely in the standard cosmological 
scenario. 

The lensing solution works better in the presence of 
a cosmological population of small dark matter clumps. 
Mini clusters of dark matter (axion-likc or Higgs-like) 
particles naturally form when a scalar field undergoes 
a second order phase transition below the QCD scale 
(Hogan & Rees 1988). Their existence does not violate 



any observational constrain provided that their mass is 
smaller than 10 4 M Q and they would behave as efficient 
lenses if their physical size is smaller than their Ein- 
stein ring, which corresponds to a lower mass limit of 
a few xl0~ 6 M Q (Zurek, Hogan & Quinn 2006). As- 
suming that all the dark matter is in mini-clusters, we 
expect that nearly 30% of the Swift afterglows could be 
strongly microlensed. To explain the observed abun- 
dance of Mgll absorbers, however, QSOs need not to suf- 
fer from these microlensing events. In other words, QSO 
beams have to be larger than the Einstein ring of the 
mini clusters (and, thus, of the GRB beams). Current es- 
timates of the QSO size (< 10 15 ~ 16 cm) then favour mini 
cluster masses which are smaller than 1 M Q . Patchy ob- 
scuration of the GRB optical beam by dust in the circum- 
burst environment could make GRB beams smaller than 
expected. Albeit very speculative, this scenario would 
certainly favour the detection of afterglows spectra with 
an increased number of absorbers even though it is dif- 
ficult to estimate the net effect on the observed dN/dz. 
Note that this result is not in contradiction with Section 
4. If the absorber size is much larger than both the beam 
sizes then the equivalent width distribution function is 
expected to be the same for QSOs and GRBs. 

An appealing possibility is that GRB afterglows are 
strongly affected by magnification bias. Light sources 
which are particularly bright in more than one waveband 
are especially likely to be lensed (Borgeest, von Linde 

6 Refsdal 1991; Wyithe, Winn & Rusin 2003). This 
"multiband magnification bias" could be very impor- 
tant for GRBs whose optical and gamma-ray luminosi- 
ties are found to be statistically independent (Nardini et 
al. 2006). It is easy to show that, independently of the 
lensing optical depth (and thus of the kind of lens), the 
fraction of lensed objects at a given redshift approaches 
unity when the mean faint-end slope of the luminosity 
functions in the two bands approaches -2 (Wyithe et al. 
2003). 7 For the faint-end of the gamma-ray luminos- 
ity function, the universal structured-jet model predicts 
a slope of 7 = —2 (Rossi et al. 2002; Zhang & Meszaros 
2002) while observational estimates based on number 
counts give 7 = — 1.57 ± 0.03 (Firmani et al. 2005) and 

7 = -1.7 ± 0.1 (Schaefer, Deng & Band 2001). The 
luminosity function of optical afterglows is not known. 
Nardini et al. (2006) find that most of the observed af- 
terglows have a similar luminosity but they do not dis- 
cuss selection effects. The fact that many afterglows are 
not detected in the optical band might suggest that the 
slope of the luminosity function is rather steep indeed. 
Magnification bias could then be responsible of the differ- 
ent counts of Mgll absorbers between GRBs and QSOs. 
Note that optically selected QSOs have a faint-end slope 
of 7 ~ -1.6 (Boyle et al. 2000; Croom et al. 2004) while 
(single-band) magnification bias becomes important for 
7 ~ —3, thereby under this scenario we do not expect 

7 The detection of an X-ray signal is often required to accurately 
locate GRB afterglows for optical follow up. Three-band magni- 
fication bias is even more efficient than the case discussed in the 
main text. For independent luminosities in the different bands, the 
critical mean faint-end slope is -5/3 (close to the observed value for 
the prompt gamma-ray emission). Note, however, that gamma and 
X luminosities of GRB afterglows seem to correlate rather tighly 
(e.g. Nardini et al. 2006) and this would again imply a critical 
slope closer to -2. 
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QSOs to be affected by microlensing events. 



5. ASSOCIATION OF Mgll SYSTEMS TO ORBS 

We now consider the possibility that the strong 
Mgll systems are associated with the GRB event. This 
solution is particularly attractive since it gives an "ad- 
ditive" correction instead of a "multiplicative" one. In 
other words, at variance with the solutions explored so 
far that multiply the incidence rate of Mgll lines by a 
constant factor, it is sufficient that a few of the 14 P06 
absorbers are associated to the GRB events to fully solve 
the discrepancy in dN/dz. Both the hint for a partial 
covering of Mgll systems along some GRB lines of sight 
and the failed optical detection of the galaxies responsible 
for Mgll absorptions (Ellison et al. 2006) support this 
hypothesis. The major limitation, however, is the need 
for cold, metal enriched gas moving at semi-relativistic 
speeds. Even assuming that the five P06 systems with 
the largest inferred velocities are produced by interven- 
ing galaxies (as expected from Sj2]), while the rest are 
associated with GRB events, we still find that the mean 
peculiar velocity of the intrinsic systems is 6 x 10 km 
s" 1 (with a dispersion of 3 x 10 4 km s _1 ). 

In the fireball model, afterglow emission is produced 
when the fastest ejecta from the central engine sweep 
up and shock the circumstellar medium. This happens 
at typical distances of 10 16-17 cm from the central en- 
gine. Gas responsible of Mgll absorption must thus lie 
further away than this. In order to produce absorption 
with column densities of 10 15 cm -2 which fully cover 
a QSO beam one needs as little as 10 -7 M© of mate- 
rial with solar metallicity. Wolf-Rayet winds, binary in- 
teractions and a SNa explosion could have easily trans- 
ported such an amount of metal enriched gas at these dis- 
tances. The GRB shock has enough energy to accelerate 
this gas to semi-relativistic velocities only within ~ 10 18 
cm from the central engine. However, most likely, the 
swept up material will also be heated to temperatures 
well above the ionization threshold of Mgll . 8 Therefore 
no Mgll absorption will be possible until the gas temper- 
ature cools down below ~ a few x 10 4 K. Even though 
the action of Rayleigh- Taylor instabilities will help re- 
shaping the structure of the clouds and produce a multi- 
phase medium, it seems unlikely that this sequence of 
phenomena could produce the observed absorption. It is 
also difficult to conceive how Mgll systems could appear 
in the intense photoionization field of the GRB afterglow. 
The characteristic isotropic equivalent luminosity of the 
optical afterglow in the Cousins R-band after 12 hours 
(in the rest frame) is 4.47 x 10 30 erg s _1 Hz -1 (Nardini et 
al. 2006). Assuming a power-law spectral energy distri- 
bution with index (3 and a time decay with index a, the 
number of ionizing photons (with energy above the sec- 
ond photoionization threshold for Mg, E > £ t hr = 15.04 
eV) emitted per unit solid angle from t mul to t = t max 
(both in seconds) is 



/43200\" 

\ ^max / 



i - 



l-a\ 



(5) 



dN> 



phot 



dn 



a 



1.82eV\ /3 

-Ethr / 



8 Unless the relativistic shock does not penetrate the dense, 
metal-rich cloudlet, which is also possible (M. Vietri personal com- 
munication) and would make Mgll absorption easier. 



For typical values a = (3 = 1.3, t min — 3 days, t max = 1 
month, this gives 10 58 photons per steradian that can 
potentially produce a second ionization of Mg atoms. 
Most of these photons will be actually absorbed by atoms 
of H,C,N and O which (assuming solar abundance ra- 
tios) are more abundant than Mg and have lower photo- 
ionization thresholds. An optically thin shell, n < 20 
nuclei per cm 3 , exposed to such a photo-ionizing flux 
gets fully ionized on extremely short timescales. Detailed 
calculations for higher density clouds require an accurate 
treatment of radiative transfer and are beyond the scope 
of this paper. The presence of structured jets (with fast 
and slow components) , already invoked to explain C iv , 
Si iv and H i absorption associated with GRB 021004 
(Starling et al. 2005), is probably needed. 

In the less favoured supranova model (Vietri & Stella 
1998), a supernova event preceeds the GRB with a time 
lag ranging from weeks to a year. Radio observations 
suggest that the fastest ejecta of "normal" supernovae 
Ibc (i.e. fast electrons) have typical velocities of ~ 0.3c 
(Berger et al. 2003) which are in agreement with the 
velocities of the absorbers in P06. If the ejecta travel 
at a constant velocity of ~ 10 5 km s^ 1 , they will need 
~ 10—100 days to cover a distance of 10 16-17 cm. Longer 
delays are needed to transport the metals produced by 
the supernova at these distances. The observed narrow- 
ness of Mgll lines (Av ~ a few hundred km s _1 ) poses 
some strong constraints on the physical size of the ab- 
sorbers. To first approximation, in an explosive event, 
Av/v = AR/R (with v velocity of the ejecta, R their 
distance from the supernova) which implies a linear size 
for the absorbers of ~ 10 12-13 cm. Thereby, these sys- 
tems could either be associated with narrow and dense 
shell-like structures or to small condensations. However, 
it is difficult to imagine how the gas could keep cold and 
to reconcile the required delays between the supernova 
and the GRB phases with observations (e.g. Meszaros 
2006). 

An alternative mechanism for producing the observed 
Mgll system requires that the material emitting the af- 
terglow is lumpy. Meszaros & Rees (1998) investigated 
spectral features arising from ultrarelativistic ions in the 
pre-afterglow fireball outflow. In this model, metal-rich 
inhomogeneities (blobs or filaments) are typically very 
small (10 5 cm) and dense (~ 10 18 cm -3 ), although their 
surface covering factor can be quite high. If these struc- 
tures survive up to the afterglow emission they could be 
responsible of the observed Mgll absorption since they 
are entrained in the flow that produces the burst. Some 
kind of magnetic mechanism should be, however, invoked 
to keep these blobs collimated in the afterglow phase. 

A promising scenario is to assume that the observed 
Mgll systems are associated with supernova remnants 
(SNRs) in the vicinity of the GRB. Mgll absorption 
lines from local SNRs and superbubbles look similar 
to those detected along high-z QSOs and appear to be 
made of several components with small velocity disper- 
sion (Danks 2000, Bond et al. 2001, Welsh et al. 2001). 
Imagine that a GRB line of sight crosses the shell of a 
young SNR that is expanding at a velocity of v ~ 3 x 10 4 
km s _1 . Also assume that the shell is spherical and 
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that it produces two absorption features (we postulate 
the presence of cold pockets of gas in the shell) in the 
GRB spectrum: one blueshiftcd and one redshifted by 
the amount (1 + zgrb)(1 ±v/c). Since GRBs are nor- 
mally associated with the highest-redshift absorbers in 
their spectra, one would erroneously assign a redshift 
(1 + zgrb)(1 + v/c) to the GRB and say that the sec- 
ond absorber moves with a relative velocity of 2v with 
respect to it. How likely is it that a GRB lines of sight 
intersects a SNR? A typical star forming region like 30 
Doradus has a characteristic size of <~ 100 pc, while a 
young SNR has a radius of ~ 10 pc. The optical depth 
for crossing a shell along a random line of sight towards 
the GRB is thus of order unity if the number density of 
SNRs is ~ 3 x 10~ 5 pc~ 3 , which corresponds to nearly 15 
objects within a single star forming region. X-ray studies 
of 30 Doradus have detected several SNR candidates and 
five superbubbles (Townsley et al. 2006 and references 
therein) thus showing that the probability for a line of 
sight to cross a SNR is not negligible. Note that not 
in all GRB spectra the expanding shell of a SNR will 
produce two absorption lines. Sometimes features of the 
host galaxy will be present together with one line from 
the SNR. In this case the velocity difference between the 
lines will be smaller. 

6. SUMMARY 

We addressed possible explanations for the different in- 
cidence rate of strong Mgll systems along GRBs which 
is about ~ 4 times higher than along QSOs. There are a 
number of phenomena that potentially contribute: dust- 
obscuration bias, different beams sizes between GRBs 
and QSOs, gravitational lensing, association of GRB ab- 
sorbers with the circumburst environment. Our results 
can be summarized as follows. 

• Considering a simplified model for Mgll absorbtion 
based on a set of realistic numerical simulations 
of galaxy formation and able to reproduce all the 
observational constraints, we showed that the inci- 
dence rate of Mgll systems in QSO spectra can be 
underestimated by a factor 1.3-2 due to dust ob- 
scuration bias. This is not enough to fully explain 
the discrepancy with the number of absorbers along 
GRBs. 

• We critically discussed the solution proposed by 
F06 where QSOs are assumed to be a factor of two 
larger than GRBs. Accounting for the dependence 
of both the equivalent width and the number den- 
sity of absorbers with the beam size, we showed 
that it is not possible to fully explain the difference 
in the observed counts by assuming that the two 
classes of background sources have different char- 
acteristic sizes. We also found that the equivalent- 
width distribution of Mgll systems detected along 
GRB lines of sight is compatible with being a ran- 
dom sampling of the QSO one at the 95.7 per cent 
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confidence level. This suggests that the absorbers 
are larger than the sources and provides evidence 
against the solution proposed by F06. 

• We showed that GRB afterglows with more than 
one absorber are brighter than the others by a 
factor of 1.7. If confirmed by increased statis- 
tics, this finding would suggest a lensing origin 
of the Mgll discrepancy. However, in the stan- 
dard cosmological scenario, lensing optical depths 
are small and difficult to reconcile with the ob- 
served counts of Mgll absorbers. This may hint 
towards the existence of an exotic population of 
lenses. For instance, microlensing due to mini dark 
matter clumps could explain the difference but this 
also requires QSO beams to be larger than GRB 
afterglows which is in disagreement with current 
observational estimates. 

• Due to their multiband selection and their emission 
properties, GRB afterglows are particularly sensi- 
tive to magnification bias. Independently of the 
nature of the underlying lenses, a large fraction of 
the observed GRBs would be lensed if the average 
faint-end slope of their gamma and optical lumi- 
nosity functions approaches -5/3 - -2. According 
to current models for prompt and afterglow emis- 
sion this is not unlikely. Present observations sug- 
gest a slope of 7 ~ —1.6 for the prompt gamma-ray 
emission-l- while the luminosity functions of optical 
afterglows is totally unconstrained. Magnification 
bias then provides a viable solution to the observed 
discrepancy. 

• The production of Mgll absorbers in the cir- 
cumburst environment (either by cold and dense 
cloudlets of metal-rich material accelerated to semi- 
relativistic speeds or, most likely, by SNRs unre- 
lated to the GRB event itself but lying in the same 
star forming region) could further reduce the sta- 
tistical significance of the observed discrepancy to 
a level compatible with statistical fluctuations due 
to small number statistics. 

We conclude that, with the possible exception of mag- 
nification bias, it is unlikely that one of these effects on its 
own can fully account for the observed counts. However, 
the combined action of some of them can substantially 
reduce the statistical significance of the discrepancy be- 
tween the incidence of strong Mgll systems in QSO and 
GRB spectra. 

We thank B. Carswell, M. Haehnelt, F. Miniati, M. 
Pettini, E. Pian, M. Rees and M. Vietri for useful dis- 
cussions. P. Madau is warmly thanked for comments on 
an early version of this paper. MV thanks the hospitality 
of ETH where this project was started. 
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